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ABSTRACT 

We study the effects of dust in hydrodynamic simulations of spiral galaxies when different 
radial metallicity gradients are assumed. SUNRISE, a Monte-Carlo radiative-transfer code, is 
used to make detailed calculations of the internal extinction of disk galaxies caused by their 
dust content. 

SUNRISE is used on eight different Smooth Particle Hydrodynamics (SPH) simulations of 
isolated spiral galaxies. These galaxies vary mainly in mass and hence luminosity, spanning a 
range in luminosities from —16 to —22 magnitudes in the B band. We focus on the attenuation 
in different wavelength bands as a function of the disk inclination and the luminosity of the 
models, and compare this to observations. 

Observations suggest different metallicity gradients for galaxies of different luminosities. 
These metallicity gradients were explored in our different models, finding that the resulting 
dust attenuation matches observations for edge-on galaxies, but do not show a linear behaviour 
in log axis ratio as some observations have suggested. A quadratic law describing the depen- 
dence of attenuation on inclination, as proposed by more recent observations, reconciles the 
attenuation of the simulations at intermediate inclinations with observations. We also com- 
pare the total infrared-to-ultraviolet flux ratios for the simulated galaxies with those of the 
SINGS sample and find general agreement. Finally we compare our results with those from 
simpler models that do not take into account structure such as spiral arms, finding that the 
inclusion of sub-structure on the size scale of spiral arms does not change conclusions about 
the attenuation dependence on inclination or wavelength. 
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1 INTRODUCTION 

Studies of dust extinction in spiral galaxies play an important role 
in this time of precision cosmology. It is remarkable how, despite 
the fact that distances have always been one of the most challeng- 
ing measurements in astronomy, we are now able to obtain accurate 
cosmological distances to galaxies and galaxy clusters. This is pos- 
sible because of relations like the Tully-Fisher relation (hereafter 
TP) and the Fundamental Plane, which relate observable parame- 
ters of galaxies, such as Hi velocity profile widths and central ve- 
locity dispersions, to their luminosity and radius. Given that these 
methods are independent of redshift, they are also important in de- 
termining the Hubble constant (Freedman et al.. 2001 ) and peculiar 
velocities of galaxies (i.e., the difference between Hubble veloci- 
ties and redshift velocities). With this, one can also obtain a poten- 
tial (.Bertschinger et al..,199Q) that can be used to compare the ob- 
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served ga laxy distributions and kinematics to cosmological predic- 
tions (e.g. lDekeletal.ll999h . Nevertheless, these valuable relations 
are subject to extinction corrections that are important to take into 
account. The scatter in the TF relation is increased if the extinction 
correction that depends on the disk inclination is not well under- 
stood. Even more so, the slope of the TF relation may be affected 
by corrections for the luminosity dependence of extinction. The de- 
pendencies of the TF relation on such corrections are demonstrated 
in detail in Giovanelli et al. (1995, hereafter G95). 

Furthermore, hierarchical models of galaxy formation predict 
that the disks of spiral galaxies form from the inside out, having 
younger stellar populations and lower metallicity in the outer disk. 
Thus, quantifying the gradients in the metallicity distribution of 
galaxies with different luminosities is important in order to com- 
pare the predictions of hierarchical models with observations. This 
can be done by looking at the effects that different metallicity gra- 
dients have on the attenuation of disk galaxies, and comparing the 
results with observations. 

Another motivation for studying the attenuation of dust in 
the disk of spiral galaxies is to gain insight into how the dust is 
distributed. The attenuation vs. inclination relations of a galaxy 
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completely free of dust and a completely opaque galaxy in which 
the obscuring dust is confined to a thin layer in the mid-plane 
of the disk are difficult to distinguish, except when seen close 
to edge-on. This is one of the reasons why determining the de- 
gree of dust attenuation in a galaxy is difficult, and has led to 
a long-standing debate a bout the extent to which gal axy disks 
are transparent or opaque I Disney et al.lll989l;[ Burstein et al]|l99ll: 
Bosma et al. 1992:'Disne v et alJl992 ':'Valentiin 1994; Kvlafis et al. 
'2001,; Holwe r da et al. lOoiT Further studies ( Moriondo et al, 1998 ; 
Xilouris et al"Hl999l) indicate that central regions of the disk may 



be opaque while outer regions are transparent. Observations re- 
veal that th e surface intensity in spirals falls off exponentially 
with radius ( de Vaucouleurs 1948) while in the vertical direction 
both the stars dvan der Kruit & Searle 1981.) and the gas and dust 
dCowie & Songailalll986l ; lKvlafis & Bahcalll 19871) fall off roughly 
exponentially with height. Exponential models of galaxies with 
equal dust and stellar s calelengths have p re viously been used to 
study dust attenuation ( Disnev et al] |l989l) . Ixilouris et alj ( [199^ 
hereafter X99) fit exponential models to observations of seven spi- 
ral galaxies seen very close to edge-on (NGC 4013, IC 2531, UGC 
1082, NGC 5529, NGC 5907, UGC 2048, and NGC 891), and re- 
vealed that in general these large spirals possess a dust scalelength 
that is about 1.4 times larger than that of the stars, and a vertical 
dust scaleheight that is about half that of the stars. These data indi- 
rectly suggest a metallicity gradient that will have different effects 
on the shape and amplitude of the attenuation vs. inclination rela- 
tion. Consequently, exploring this relation can lead to conclusions 
about the radial distribution properties of the dust in the disk of 
spiral galaxies. 

Other observations have quanti fied the metall i city g radient of 
spiral galaxies by different methods. IZaritskv et al. I ( 1 19941 . hereafter 
Z94) provides the metallicity gradient of 39 disk galaxies, by exam- 
ining the properties of HII regions within the disk of these galaxies. 
They found an average metallicity gradient of —0.05 dex/kpc, al- 
though this value varied widely from 0.00 to —0.20 dex/kpc. They 
found a correlation between abundance gradients and Hubble type 
only when the gradients are expressed in terms of dex/kpc as op- 
posed to dex/isophotal radius. 

Based on the fact that colour gradients in disk galaxie s are due 
in most part to the gradients in age and metallicity ( Bell & de Jon^ 
l2OO0h . lMac Arthur et al] ( l2004 hereafter MA04) provided results of 
metallicity gradients from observations of optical and NIR galaxy 
colour gradients in a sample of 172 galaxies with different lumi- 
nosities. They found that the metallicity gradients in their sample 
were steeper than those of Z94 for galaxies fainter than ~ —20, 
and shallower for galaxies brighter than ~ —20, in the K-band. 

G95 sampled more than 1700 galaxies for which I-band CCD 
images, redshifts, and HI line widths were available. Their conclu- 
sions were that their attenuation as a function of inclination could 
be well fitted by a linear function of log(a/6), where a/b is the 
axis ratio. From sampli ng 87 spiral galax ies within the Ursa Ma- 
jor and Pisces Clusters, ITuIIv et aU ( Il998l . hereafter T98) reached 
conclusions si milar to those of G9 5, with their data fitting well a 
linear relation. [Masters et alj ( |2003| . hereafter M03) based their re- 
sults on data take n from 3035 galaxies fro m the 2 Micron All-Sky 
Survey (2MASS Ijarrett et alJboOll , l2003h for which I-band pho- 
tometry and redshifts were available, and concluded that a simple 
linear law was not adequate for describing the attenuation depen- 
dence on inclination observed in their sample. Thus, M03 adopted a 
bilinear law fitting separately high and low inclinations systems as 
separated by log(a/6) = 0.5. Looking for even less linear shapes 
that fit their observed colour dependence on inclination, M03 also 



tried a quadratic law. The results provided by M03 and G95 cover 
only the infrared and near infrared range of the spectrum. 

Here we use the capabilities of SUNRISE, a Monte Carlo 
radiative-transfer code t Jonsson 2006), to perform "simulated ob- 
servations" of dust attenuation in the disks of eight different hydro- 
dynamic simulations of spiral galaxies, having different luminosi- 
ties and sizes, while assuming the different metallicity gradients 
provided by MA04, X99, and Z94. Performing these "simulated 
observations" at different inclinations for each of the galaxy mod- 
els, we obtain the attenuation dependencies on inclination and lu- 
minosity for each of the assumed metallicity distributions. These 
results are then compared to the observations made by M03, G95, 
and T98. By finding which simulations best fit the data, it is de- 
termined which metallicity gradients are the most appropriate for 
these galaxy models. 

The main difference between this work and previous ways of 
quantifying the attenuation produced by dust in galaxies through 
radiative transfer calculations is that our calculations are per- 
formed over an arbitrary geometry specified by the hydrody- 
namic simulations. Thus our galactic models resolve structure such 
as spiral arms that in previous radiative transfer models, either 
by analytical appro ximations (^ yun et al. 1994; Silva et al. 199^; 
Ixilouris et al][l999l ; Isaes & Dei onghe 2001; T uffs et alJi2004h or 
by Monte Carlo c alculations Ide Jong,. 1996 ; Kuchinski et a l .11 19981; 
Ferraraet alj[l99 9; Bianch i et alJI 19961 [2OOO; Matthews & Woo j 



2001l ; lPierini et aU200 4). is not considered. Nevertheless, the struc- 



ture resolved in our simulations is not small enough to include in- 
dividual star forming regions, which affects the resulting attenua- 
tion. We compare our results with previous models to examine if 
the inclusion of structures such as spiral arms makes a difference 
in the dependence of attenuation on inclination and wavelength. Fi- 
nally we have extended the triple exponential model o f Disnev et al] 
( Il989l) to explore how the dependence of attenuation on inclination 
changes with different dust-to-star scalelengths ratios. 

The structure of this paper is as follows: Section 2 is a presen- 
tation of the hydrodynamic galaxy models. In section 3, a descrip- 
tion is given of how the simulated observations were made with 
SUNRISE. In section 4 the results are compared to real observa- 
tions. In section 5 we compare to previous models of galactic disks 
without structure such as spiral arms, and explore how the depen- 
dence of attenuation on inclination changes with different dust-to- 
star scalelength ratios. Section 6 provides a summary and lists the 
conclusions. 



2 THE HYDRODYNAMIC GALAXY MODELS 

Radiative-transfer calculations were made for eight different 
Smooth Particle Hydrodynamics (SPH) simulations of isolated 
galaxies. These simulated galaxies were created using a ver- 
sion of the code GADGET 1 which enables entr opy conservation 
dSpringel et alj200ll ; [Springel & Hemauisll2002h . as part of a com- 
prehensive work aimed at studying galaxy merge rs (Cox 20o3 
ICox et al ]|2006l. I2OO7I ; IJonsson et alj2006l ; |j onsso 3r2006h . 

The set of isolated galaxies for which calculations of dust ex- 
tinction were made includes four models that approximate late-type 
spiral galaxies and four less gas-rich spiral galaxies (the G-series). 
The late-type spirals group comprises an She type galaxy, a more 
massive and bigger Sbc (Sbc+), a less massive and smaller Sbc 
(Sbc-), and an Sc type. The G-series have been created spanning a 
wide range of masses and having average properties for their mass. 
Table [T] lists the characteristic of each model. The details of these 
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!;alaxy models are described in ICox etaljJioO^ and ljonsson et alj 
20061) . Here we mention some of their important features. 

The galaxies contain a stellar and gaseous disk, which are 
rotationally supported, a spheroidal stellar bulge, and a massive 
dark matter halo . A model for supernova feedback, as described in 
ICox et aL I l l2006h , is included and is crucial for the stability of the 
gas disk and for preventing excessive star formation. The adopted 
star f ormation prescription is based on a Schmidt law ( Schmidt) 
Il959[) and is consistent with the Kennicutt law ( Kennicutt 1998). 
Schemes for radiative cooling and metal enrichment have also been 
implemented in these models. Our models do lack age gradients, 
something that would be necessary for comparing our results with 
the predictions of hierarchical models of galaxy formation. 



3 THE RADIATIVE TRANSFER CALCULATIONS 

3.1 SUNRISE 

SUNRISE is a high-performance code designed to combine a full 
radiative-transfer model with high-resolution hydr odynamic simu- 
lations of either isolated or interacting galaxies; see IjonssonI f20o3) 
for a detailed description of how SUNRISE achieves this goal. SUN- 
RISE first assigns an initial spectral energy distribution (SED) for 
the stellar particles in the hydrodynamic simulations. Assuming 
that the density of dust traces the density of metals, it then solves 
the radiative-transfer problem by the Monte-Carlo method in the 
given arbitrary geometry for a number of wavelengths, and fi- 
nally, it interpolates to obtain a resulting full SED after absorp- 
tion and scattering. The SEDs used are from Starbust9^Q v5.0 
jLeitherer et al.|[l999) with Kroupa IMF. In post-processing, SUN- 
RISE calculates total magnitudes and generates images for the spec- 
ified optical filter bands, as well as the total absorbed bolometric lu- 
minosity, over all wavelengths. In this way SUNRISE can be used to 
make simulated observations that can be easily compared to actual 
observations. SUNRISE is free softwar^f]. 

For this investigation, calculations have been done using eight 
different filters (U, B, V, R, I, J, H and K) and thirteen cameras at 
different inclinations. Six cameras were placed above the plane of 
the disk at inclination angles 0, 35, 51, 63, 74, and 85 degrees. 
Another six were placed at the same inclinations but on the op- 
posite side of the disk and at a 90 degree azimuthal angle to the 
first set. One more camera was placed exactly edge-on. These cal- 
culations were done for two simulation times, 0.5 Gyr and 1 Gyr. 
This allows enough time for the simulation to relax but is not so 
late that a significant amount of gas is consumed. For each inclina- 
tion, the results were averaged over the two azimuthal angles and 
the two snapshots. For simplicity, the dust model used in all our 
models regardless of s i ze wa s the R = 3.1 Milky- Way model by 
IWeingartner & Draind fcoOlh . 

3.2 Star-Formation History and Metallicity Distributions 

Given that the metallicity and the age are required to pick the SED 
of a stellar particle, two important values that are given as inputs to 
SUNRISE are the star-formation history and the metallicity distribu- 
tion of the stars present when the simulation is started. 

In all the calculations made for this paper the bulge stars are 
assumed to have formed in an instantaneous burst, while the disk 

|http : / /www ■ stsci ■ edu/science/starburst99/| 
http : / /sunrise, f amil jen jonsson . org| 



has had an exponentially declining star-formation rate starting at 
the time of bulge formation and leading up to the start of the simu- 
lation. The age of the bulge and stars at the beginning of the simu- 
lation, and the time constant r of the exponential decline used for 
each of the galaxy models are given in Table [T] The time constant 
r for each galaxy was determined by fitting an exponential star 
formation rate history to the following constraints: the present-day 
star-formation rate must match the star formation rate in the hydro- 
dynamic simulations, the total stellar mass must match that in the 
simulated models, the average stellar age must match the data in 
MA04 (as well as possible), and finally no stars can be older than a 
Hubble time. In the case of the Sc galaxy, this means a marginally 
growing star-formation rate (at the onset of star formation 13.8 Gyr 
ago, the star formation rate was 37% of the present-day SFR). This 
is not very surprising, since a large, gas rich spiral today cannot 
have sustained its large star formation rate for a long time without 
exceeding its current stellar mass, it's just being assembled later 
than average. Effectively, the very large numbers (positive as well 
as negative) just mean a constant star formation rate history. 

The metallicity of the gas and stars at the start of the simula- 
tion is, like the surface densities, assumed to decline exponentially 
with radius in the disk. The metallicity of the gas particles is as- 
sumed to be the same as that of the stars and thus the same depen- 
dence applies for the gas metallicity as a function of radius. lOwekl 
(1998) concluded that the formation of dust grains in the interstellar 
medium is, in good approximation, linearly proportional to the den- 
sity of metals, with a proportionality constant of 0.4. If that is the 
case, the dust density d{r) would also decline exponentially from 
the centre as is observed by X99, and be described by a relation of 
the form 

d{r) = OAg{r)z{r) = doe-'^-^nO"^" = d^e-"/^-* , (1) 

with Rg and Rd being the gas and dust scalelengths respectively, 
do — OAzogo being the central density of dust, and G being a con- 
stant that determines the gradient of a metallicity distribution ex- 
pressed in a decimal exponential manner, as is usually given (Z94, 
MA04). The observations of X99, as described in Section 1, sug- 
gest that the dust scalelength of big spiral galaxies is about 1.4 
times greater than the stellar scalelength. Hence, we can define the 
"Xilouris gradient" Gx as the gradient that results in a ratio of dust 
scalelength to stellar scalelength equal to 1.4. Recalling equation[T] 
we have that 

or 

The Xilouris gradient Gx for each of the galaxy models is shown 
in Table |2] 

X99 observed only big spiral galaxies, so Gx is not necessar- 
ily the best metallicity gradient assumption for our set of smaller 
galaxies. MA04 (Fig. 26) provides metallicity gradients as a func- 
tion of luminosity for a sample of 127 galaxies, from which we 
have obtained another set of gradients that were applied to the three 
smallest galaxies, G2, Gl and GO. Table [2] shows these last gradi- 
ents and summarises all the gradients that were tested in the simu- 
lations. 

It is important to note that the central metallicity of the galax- 
ies depends on the gradient in use. This is because the metallicity 
given in Table[T]is at 1.3 scalelengths from the centre. 
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Table 1. The galaxy models used for the calculations of dust attenuation. 



Model M^ir^' Afft^ ij^c 2^//?/ Rg/Rf /g* R^'^ y^t' ^1.3^ Ng^ Age' r 

(M0) (Mq) (kpc) (kpc) km/s (Zq) (xlO*) (Gyr) (Gyr) 



Sbc+ 


9.28 ■ 


10" 


1.56 ■ 


10" 


7.0 


0.125 


3.0 


0.52 


0.10 


0.60 


210 


1.12 


3 


13.9 


110 


Sc 


8.90 ■ 


10" 


1.12 ■ 


10" 


4.7 


0.2 


4.0 


0.69 


0.00 


0.00 


196 


1.00 


3 


13.8 


-13 


She 


8.12 ■ 


10" 


1.03 ■ 


10" 


5.5 


0.125 


3.0 


0.52 


0.10 


0.45 


195 


1.00 


3 


13.9 


-106 


G3 


1.16 ■ 


10^2 


6.22 ■ 


IQio 


2.8 


0.125 


3.0 


0.20 


0.14 


0.37 


192 


1.00 


5 


14.0 


10 


Sbc- 


3.60 ■ 


10" 


4.98 ■ 


IQio 


4.0 


0.125 


3.0 


0.52 


0.10 


0.40 


155 


0.70 


3 


13.7 


124 


G2 


5.10 ■ 


10" 


1.98 ■ 


IQio 


1.9 


0.2 


3.0 


0.23 


0.08 


0.26 


139 


0.56 


3 


14.0 


8.2 


Gl 


2.00 ■ 


10" 


7.00 


■ 109 


1.5 


0.2 


3.0 


0.29 


0.04 


0.20 


103 


0.40 


2 


11.5 


3.7 


GO 


5.10 ■ 


lOio 


1.60 


■ 109 


1.1 


0.2 


3.0 


0.38 


0.01 


0.15 


67 


0.28 


1 


8.7 


1.4 



''Virial mass. 
''Baryonic mass. 
'^Stellar disk scalelength. 

Ratio of stellar-disk scaleheight and scalelength. 
° Ratio of scalelengths of gas and stellar disks. 
'Gas fraction (of baryonic mass). 
^Bulge fraction (of baryonic mass). 
'^Bulge scale radius. 
'Circular velocity. 

j Metallicity at 1.3 scalelengths from the centre (gas and stars), from Z94. 
'^Number of gas particles. 

'Age of oldest stars (formation time of bulge and oldest disk stars). 
"^Exponential time constant of the star formation rate for the disk stars. 



Table 2. Tested Metallicity Gradients. 



Galaxy 


Z94 


X99 (Gx) 


MA04 




{—dex kpc~^) 


{—dex kpc~^) 


{—dex kpc~^) 


Sbc-l- 


0.05 


0.023 




Sc 


0.05 


0.043 




Sbc 


0.05 


0.030 




G3 


0.05 


0.058 




Sbc- 


0.05 


0.041 




G2 


0.05 


0.086 


0.04 


Gl 


0.05 


0.112 


0.05 


GO 


0.05 


0.148 


0.06 



The adopted metallicity gradients in our models are shown in bold. 



The advantage of producing simulated observations with SUN- 
RISE is that the results of the calculations can be directly compared 
to observations, which will be done in the next section. 



4 RESULTS AND COMPARISON WITH OBSERVATIONS 

Before analysing any data, we can anticipate some of the results 
by simply looking at the colour composite images obtained in our 
calculations. Figure [T] clearly shows the effects of having or not 
having a metallicity gradient when the Sbc model is seen edge-on, 
while Figure|2]shows how the metallicity, and hence the strength of 
the dust lane an d atten uation, decrease with luminosity. Note that 
iDalcanton et al.l ( |2004|) found that galaxies with rapidly rotating 
disks (i.e. Vrot > 120 km/s) have well defined dust lanes, while 
galaxies with lower rotation velocities do not show evident dust 
lanes. Figure |2] shows that it is hard to see any dust lanes in our G2 
and smaller galaxy models, which have rotation velocities of less 
than 139 km/s (See Table[T]l. 



4.1 Observations of Disk Attenuation. 

Some of the results of G95, T98, and M03 will be used to compare 
the data obtained from our simulations with observations. These 
authors have quantified the extinction effects in the disks of galax- 
ies as a function of inclination, with an aim to better understand the 
necessary extinction corrections to the TF relation. The extinction 
effects of dust are normally quantified by measuring the attenua- 
tion, which refers to the decrease in flux resulting from photons 
being scattered and/or absorbed. This is defined as 

A>. = ml -mx^ -2.5 log (^^^ , (4) 

with m\ and /° being the intrinsic magnitude and flux in the ab- 
sence of dust respectively. In the simulations, it is simple to find the 
attenuation dependence on inclination of a galaxy by simply com- 
puting the attenuation, as in equation ID for a number of different 
lines of sights at different inclinations. In real systems, it is not easy 
to determine the actual inclination of the disk and it is impossible to 
know the intrinsic magnitude or flux by direct measurement, hence 
more practical ways of defining inclination and attenuation are re- 
quired. 

The technique used (G95, T98 and M03) to find the inclina- 
tion is by first determining an isophotal radius (a radius of constant 
surface brightness) in the disk of the observed galaxies, and then 
finding the ratio between the semi-major and the semi-minor axis 
of the isophotal ellipse, which is related to the actual inclination by 

cost = J , (5) 

where a and b are the semi-major and semi-minor axis respectively, 
and q is the axial ratio when the system is viewed edge-on. Having 
this and given that surface brightness is proportional to the flux, one 
could expect that, as a simple approximation, the attenuation rela- 
tive to face-on as a function of axial ratio would follow a relation 
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Figure 1. Colour composite images (I,V,B bands) of tlie Sbc model with 
no dust, assuming the Xilouris gradient, and assuming constant metallicity 
for top, middle, and bottom respectively. We can see a significant difference 
between these three cases when edge-on. 

of the form 

^ J / face-on area = na^ \ ^ / a \ 

A 7a g ^ inclined area = nab J "^^ ^ \ 6 / ' 

where 7^ is an attenuation parameter as a function of wavelength 
A and luminosity. Of course this would be completely valid only if 
the surface brightness of the system were constant. G95, T98 and 
M03 fitted their observations with a relation of the form of equation 
|6]and with that obtained the most appropriate value of 7^. 

Since it is not possible to know the intrinsic magnitudes and 
these observers were only interested in the dependence on inclina- 
tion of the attenuation in a given band, the observations used for 
their fits were those of the A — K colour (where A can be any given 
band). This is because dust is optically thin in the near infrared part 
of the spectrum, and thus the magnitude in K should be approxi- 
mately constant with inclination. Thus, the main technique used for 
these observations is to find the 7a that best fits their observations 
of the dependence of the colour A — K on the axial ratio a/b. Under 
the assumption that extinction effects in the K-band are negligible 
and independent of inclination, the dependence of A — K colour on 
inclination would then represent the actual attenuation dependence 
on inclination for the given band A. 

T98 assumed that the extinction effects in the K-band, al- 
though small, were not entirely negligible. They concluded that 
their results of attenuation as a function of inclination were well 
fit by a linear relation in logarithmic space that followed the form 
of equation [6] (see Fig 2 in T98), with a slope 7a that is dependent 
on luminosity. Their adopted equations for 73 and 71 were 

7b = -0.35(15.6 + Mb +5 log /180) (7) 



Figure 2. Face-on and edge-on colour composite images (I,V,B bands) of 
the Sc, G3, G2 and Gl galaxy models from top to bottom respectively. The 
scale is the same so that the differences in size can be visualised. The GO, 
which is not shown here, is even smaller and less luminous than the Gl 
(see Table [T). Notice how the strength of the dust lane and the reddening 
decreases with the size (luminosity) of the galaxy, and is unnoticeable for 
the Gl and G2. 

71 = -0.20(16.9 + Mi +5 log /i8o) , (8) 

with Mb and Mi being the magnitudes corrected to face-on for 
the B-band and I-band respectively, and Hq = /igoSO km/s/Mpc. 
The conclusions of G95 were the same as those of T98, with a fit 
described by equation[6l and a luminosity dependence given by 

71 = -0.16(15.25 + A/i). (9) 

However, M03 concluded that a simple linear law as defined in 
equation[6]was not adequate for describing the dependence of atten- 
uation on inclination observed in their sample. Thus, M03 adopted 
a bilinear law fitting separately high and low inclination systems as 
separated by log(a/&) — 0.5. Looking for even less linear shapes 
that better fit their observed colour dependence on inclination, M03 
also tried quadratic fits of the form 

Ax.^Bx. + Ca. log(a/6) + Da. [log(a/6)]' , (10) 

where C\v, B\^, and Da. are constants that best fit their observed 
A — colour dependence on axial ratio. 

4.2 Comparison of Results with Observations 

First, the relative attenuation between the K-band and some other 
band was examined. As mentioned in the last sub-section, T98 as- 
sumed that the extinction effects on the K-band, although small, 
were not entirely negligible. They assumed that, independent of in- 
clination, Ak/Ab relative to face-on was 0.15, and Ak/Ai was 
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Table 3. Ratios of Relative to Face-on Attenuations. 
e(deg) X99 Z94 X99 Z94 


















35.4 


0.30 


0.42 


0.09 


0.20 


50.7 


0.30 


0.40 


0.14 


0.21 


63.2 


0.39 


0.35 


0.21 


0.20 


74.4 


0.34 


0.35 


0.20 


0.23 


85 


0.38 


0.39 


0.26 


0.28 


90 


0.41 


0.43 


0.30 


0.32 


Average 


0.39 


0.33 


0.20 


0.24 



The ratios of relative to face-on attenuations between the K. I and B 
bands for the Xilouris and Zaritsky gradients at different inclinations 
for the Sbc galaxy model. 



0.25. Our results for these ratios with both the X99 and Z95 gradi- 
ents are given in Table [3] 

Table [3] shows that the average relative attenuations found in 
our simulations when the Xilouris gradient is assumed (0.20 and 
0.39 for Ak/Ab and Ak/Ai respectively) are just a little higher 
than the ratios assumed by T99. However, these ratios vary signif- 
icantly and in a non-trivial manner with inclination. This was not 
considered by T98. If this fact is taken into account, their values of 
observed attenuations should be about 15% higher for large edge- 
on galaxies. For smaller galaxies and intermediate inclinations the 
difference is not as large. 

Figure [3] shows a direct comparison of the attenuation depen- 
dence on inclination and on luminosity with observations. For the 
luminous galaxies, the Sbc variants, the Sc, and the G3, the Xilouris 
gradient results in edge-on attenuations which are roughly consis- 
tent with the observations. However, this gradient is too steep for 
the low-luminosity galaxies, producing too little attenuation com- 
pared to observations, and X99 only studied large spirals. For this 
reason, the metallicity gradient for the G2, Gl and GO galaxies was 
taken from the results of MA04. In some cases, notably the G3, G2 
and G 1 , the no-gradient simulations seem to match the observations 
best, but as galaxies are observed to have metallicity gradients, this 
model is not acceptable and we chose the simple alternative of us- 
ing the X99 gradient for the luminous galaxies and the MA04 gra- 
dient for the smaller ones. 

It is apparent from Figure [3] that, although the simulations 
match some of the observed amplitudes, the dependence of atten- 
uation on inclination is certainly not linear; a relation of the form 
Ax = 7a log(a/6) is not appropriate. The bilinear law adopted by 
M03 (dashed green in Figure [3} suggests a less linear relation, but 
still is not similar to what we obtain in the simulations. We also 
compared the M03 quadratic fits with our simulations. M03 (see 
Figure 4 in M03) used the statistical average of the colour depen- 
dence on inclination observed in their full sample to obtain the fits 
shown in Figure|4l In order to compare with them, we observed the 
colour dependence on inclination for each of our eight simulated 
galaxies, using the adopted gradients for each (see Table[2}; this is 
shown in Figure|4] The shape of the quadratic fits agrees rather well 
with the results of the simulations. 

Lastly, we compared the attenuation dependence on luminos- 
ity in our galaxy models with what was observed by G95 and T98. 
T98 gives results of luminosity dependence for both the I and B 
bands (equations |7] and [8]l, while G95 only gives results for the I- 
band (equation[9]l. This is shown in Figure|5] which illustrates that 
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Figure 4. The colours J-K, J-H, and H-K relative to face-on plotted against 
the log of the axial ratio. The black dashed lines represent the quadratic fits 
described by equation [To] Over-plotted are the colour dependencies found 
in our models. We do not make a best fit to our points due to the fact that 
we have more bright galaxies than faint galaxies. The shape between simu- 
lations and observations is very similar, and a relation of the form given by 
equation [Tol accuratelv describes the attenuation dependence on inclination 
in our simulations. 



our points closely follow the fits of G95 and T98, supporting the 
conclusion that the adopted metallicity gradients are good for our 
galaxy models. These gradients and the face-on luminosities of our 
galaxy models when using these gradients are given in Tables|2]and 

s 

So far we have been examining the attenuation in certain filter 
bands. However, SUNRISE calculates dust attenuation at all wave- 
lengths, so we can also calculate the total bolometric stellar lumi- 
nosity that is absorbed by dust and consequently re-emitted in the 
infrared. When mid/far-IR observations are available, this "bolo- 
metric attenuation" provides a sensitive measure of the total amount 
of dust, thus avoiding the problem of determining the face-on atten- 
uation. 

In order to compare wit h the Spitzer Infr ared Nearby Galax- 
ies Survey (SINGS) sample jPale et aljr2007l hereafter D07), we 
have calculated the infrared-to-ultraviolet ratios in our models in 
the same manner as they did, by dividing the total infrared lumi- 
nosity (TIR) with the ultraviolet luminosity given by the sum of the 
luminosities at 1500 A and 2300 A (i.e. the FUV and NUV filter 
bands of the GALEX satellite). Figure |6] shows the results for our 



Dust Attenuation in Spiral Galaxies 7 



B-band 



l-band 



2.0 

1.5 

1.0 

0.5 

0.0 
2.0 

1.5 

1.0 

0.5 

0.0 
2.0 

1.5 

1.0 

0.5 

0.0 






■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ] I ' ; 

no-gradient / " 

X99 / : 
Z94 / 
G95 / 
T98 / / ■ 
IVI03 /-Sy-y ' 


^ — , — , — 1 — , — , — , — 1 — , — , — , — 1 — , — , — ^pH — , — ,— 

: Sc / - 


— , — , — 1 — , — , — , — 1 — , — , — , — 1 — , — , — , — h— ■ — .-^ 


— , — , — , — 1 — , — , — , — 1 — , — , — , — 1 — , — , — t-H — 1 — 1— 

: Sbc / / ■ 


— , — , — , — 1 — , — , — , — 1 — , — , — , — 1 — , — , — , — 1— , — ,— • 



0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 

iog(a/b) 



Figure 3. Plots of attenuation relative to face-on vs. the log of the axial ratio, in the B-band (left column) and the I-band (right column), for our three most 
luminous galaxy models: the Sbc+, Sc, and Sbc (top to bottom). The solid lines represent the simulations for two metallicity gradients, Z94 (green) and X99 
(blue), and a run with no gradient (red). Dashed lines are the best fits to the observations of G95, T98 and M()3. The green dashed Une is the biUnear line 
adopted by M03, valid only for their most luminous galaxies, while the rest are given by equations[2][l]and|9]varying with luminosity. Although it is not evident 
if one of the gradients matches better the amplitude of the observations in I-band, it is noticeable that the X99 gradient matches the edge-on attenuations better 
than the Z94 gradient in the B-band (especially considering that T98 observations should be about 15% higher due to K-band attenuation), hence we adopt 
it as our preferred gradient for the most luminous galaxy models. Our simulations evidently do not match the shape of the fits adopted by observers, as they 
predict less attenuation for intermediate inclinations than the observations show. 



Table 4. Luminosities 



Galaxy 


Mb 


Mi 


Mk 


Labs 1 Lhol 


Sbc-l- 


-21.6 


-22.9 


-23.2 


0.48 


Sbc 


-21.4 


-22.5 


-22.8 


0.44 


Sc 


-21.2 


-22.6 


-22.6 


0.46 


Sbc- 


-20.3 


-21.6 


-21.9 


0.37 


G3 


-20.4 


-22.0 


-22.4 


0.32 


G2 


-19.2 


-20.8 


-21.2 


0.32 


Gl 


-17.9 


-19.6 


-19.9 


0.14 


GO 


-16.0 


-17.6 


-18.1 


0.02 



Face-on luminosities in B, I, and K bands for our set of galaxy models 
when the adopted metallicity gradients of Table[2]are used. The last column 
shows the fraction of bolometric luminosity which is absorbed by dust. 



different models at different inclinations. We only compare within 
the range of inclinations included in the DO? sample (i.e. to 70 



degrees). Comparing the range of infrared- to-ultraviolet ratios for 
the Sbc, G3 (which approximates an Sb type), and Sc models with 
the galaxies of these types observed by D07 (see Figure 5 in DOT), 
the simulated galaxies show slightly lower attenuation. While DO? 
observed ratios from about 1 to 10, we obtain ratios from about 1.5 
to 7, but the SINGS sample only contains a few star-forming galax- 
ies with these morphological types. Overall, the amount of dust at- 
tenuation in our galaxies, as given by the infrared-to-ultraviolet flux 
ratio, appears to be fairly close to what is observed in real galaxies. 

Another measu r ement of the bolometric attenuation was done 
bv lPopescu & Tuffs! j2002l) using the ISOPHOT instrument on the 
ISO satellite, who found correlations between the overall attenua- 
tion and the Hubble type of a sample of galaxies in the Virgo clus- 
ter. They obtained a mean bolometric attenuation in galaxies of type 
Sab-Sc of about 0.23. Comparing the val ues in Table[4]for the G3 , 
Sbc-, Sbc, Sbc-l-, and Sc to the results of IPopescu & Tuffj ( l2002h . 
the simulations have attenuations that are about 15 — 20% higher. 
(The simulated galaxies here have slightly lower bolometric attenu- 



8 Rocha et al. 



B-band 



l-band 



2.0 

1.5 

1.0 

0.5 

0.0 
2.0 

1.5 

1.0 

0.5 

0.0 
2.0 

1.5 

1.0 

0.5 

0.0 

1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
0.4 
0.3 
0.2 
0.1 
0.0 




^^^^^^^^ 


■ 1 ■ 1 1 1 ■ 1 1 ■ ■ 1 1 ■ 1 1 ■ ■ 

: no-gradient / : 

X99 / 
Z94 / : 

' IVIA04 / " 
G95 // : 
T98 / 


• — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — - 

- G3 : 

••-''''''''''''^ 


■ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — . 
J 


• — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — - 
: G2 : 

^-^-^-=====^^^ : 


; 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ■ 






■ — \ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — - 
■ G1 j 


■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ■ 

_ 






: 1 1 1 — 1 1 1 1 — 1 1 1 1 — 1 1 1 1 — 1 1 1 ; 


\ — ' — 1 — ' — 1 — 1 — 1 — ' — 1 — 1 — ' — ' — 1 — 1 — ' — 1 — 1 — ' — ' — : 







0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.1 
log(a/b) 



Figure 3 - continued The Sbc-, G3, G2, Gl, and GO (top to bottom). We can see how, for the G-series, the relative attenuation with the X99 (solid blue) 
gradient falls below the observations. Of the observed metallicity gradients, those of MA04 (solid pink) give a relative attenuation that is closest to what 
is observed, hence we have adopted these as our prefeiTed gradients for the GO, Gl, and G2. The results using the Z94 (solid green) for the GO, Gl, and 
G2 are similar and overlap in some cases with those of MA04. For the lower-mass galaxies, using constant metallicity is sometimes the best-fitting model. 
Nevertheless, the no-gradient model was not adopted since no observations support such a metallicity distribution. 



Dust Attenuation in Spiral Galaxies 9 





' 1 






1 1 i i 1 




' 1 


1 1 1 1 1 


' 1 


1 1 i 1 i 1 i 
















T98 




GO + - 


4 














G94 




G1 )K : 




















G2 □ 




















G3 O ■ 




















Sbc A 




















Qhr- 1 — 1 
OUU 1 1 


3 


















Sbc+ X 




















Sc ♦ 




















A 


2 














































□ 


o 










o ■ 


1 






































□ ] 

























: .+^-'' 
■ 1 


1 









1 


-■■■■+ 
, -■; , 1 


1 





j: I I 1 1 I I I I I I 1 1 I I I ^ I . --. I I I 1 I I I I 1 I 1 I 1 I 

-16 -18 -20 -22 -16 -18 -20 -22 -24 



magnitude magnitude 



Figure 5. The attenuation parameter 7, as defined by the relation |6] plotted against the face-on magnitudes of all our models in the B-band (left), and the 
I-band (right). The dotted and dashed lines are the best fits to the observations of T98 and G95 respectively. Our points match the I-band observations less 
well, but the agreement is still reasonable given the fact that we only have two galaxies fainter than —20. 
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Figure 6. Infrared-to-ultraviolet ratios for our models, as a function of 
inclination. This figure can be directly compared to Figures 4 & 5 in 
D07. Compared to the SINGS galaxies of the coiTesponding morpholog- 
ical types, the simulations show slightly lower dust attenuations but cover 
roughly the same range of values. The total infrared emission in our models 
thus appears to be consistent with what is observed by SINGS. 



adons than those previously obtained bv ljonsson et ai]l l2006h usms 

simulations with constant metallicity.) 

Apparently, the results of D07 and IPopescu & Tu^ ilOO'i) 
are mutually inconsistent. It is possible that the ISOPHOT obser- 
vations, only sensitive to wavelengths longer than 60 fim, missed a 
significant fraction of the dust luminosity at shorter wavelengths. 
Another possibility is that the spiral galaxies in the Virgo clus- 
ter have systematically lower gas content and lower star-formation 
rates than the average field spiral. Overall, the fact that our galaxies 



have slightly lower bolometric attenuation compared to the newer 
and higher-quality data of the SINGS galaxies, as well as slightly 
lower values of , especially at intermediate inclinations, paints a 
consistent picture that the simulations have dust attenuations on the 
slightly low side of observed galaxies of the same morphological 
types. 



5 COMPARISON WITH OTHER MODELS 

A simple way to model dust attenuation in the disk of galaxies 
is to perform radiative transfer calculations in smooth exponen- 
tial distributions of stars and dust, disregarding scattering. This 
is unrealistic but helps to understand the behaviour of attenuation 
when changing parameter s such as in c linatio n or dust-to-stellar 
scalelengths/scaleheights. ,Disnev et al. I ( 1 19891) developed a triple 
exponential model in which a single scalelength describes the ra- 
dial distribution of stars and dust, and two vertical scaleheights, 
one for dust and one for stars, describe the vertical distributions 
of the model. This sa me model was later used and extended by 
iGiovanelli et alj ( Il994i) and M03. Motivated by the results of X99, 
we explored the same model but allowing different scalelengths 
for dust and stars (a "quadruple exponential" model). In a model 
like this, the emissivity and absorption coefficients of the radiative- 
transfer equation would be given by 



e(r,z) 



K(r, z) 



(.0 exp 



Ko exp 



r 
r 



El 

Zs 

N 



and 



(11) 



(12) 



with £0 and Ko being the emissivity and the absorption coefficient at 
the centre, respectively, Rs and Zs the scalelength and scaleheight 
of the stars, respectively, and Rd and Zd the scalelength and scale- 
height of the dust, respectively. Figure|7]plots the resulting attenua- 
tion dependence on inclination for different values of A = Rd/Rs, 
and with ^ = Zd/Zs = 0.5, q = Zs/Rs = 0.15, and t = 1, 
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Figure 7. Attenuation relative to face-on vs. the log of the axial ratio for 
the "quadruple exponential" model when different values of A = / Rs 
are given. We can see that the different values of A result in different mag- 
nitudes of edge-on attenuations, but do not affect the quadratic shape of the 
relation much. 



where r is the central face-on optical depth 



«;(0, z) Az = 2ko 



Az = 2KoZd. (13) 



iDisnev et alj h989l) , lGiovanem et alj fl994h , and M03 show how 
this simple model behaves when different values of ^, q, and r 
are given. We find that when different values of A are assumed the 
edge-on attenuation varies significantly, but the shape of the curve 
remains the same (see Figure|7](. 

More sophisticated radiative transfer models where 
anisotropic scattering is considered and different parameterisations 
and spatial distributions of clumps are included c an be found 
extens ively in the literature (s ee the introductions of IPierini et alj 
| |2004|) and iTuffs et id] ( l2004h for nice reviews of the history of 
these models), but few explore the consequences of structure 
such as spiral arms in the resulting attenuations jCorradi et alj 
1 19961 : ISemionovet"ai]|2007h . If we compare our results for the 
attenuation depend ence on inclinat i on (F ig ure [3) with s o me of 
these models, i.e. 'Giovanelli et al ('1994'); 'Bianc hi et al] jl996l 
[2OOO); Ferrara et al. ( 1999); Baes & Dejonghe (200lh; lBvun et alj 
d 19941) ; ITuffs et alj (12004) we cannot see any significant difference 
that can be attributable to the geometry of our simulated galaxies. 

The dependence of attenuation on wavelength in our Sbc 
model (which has a central optica l depth r close to 2) has been 
compared with that of lPierini et alj J2004. hereafter P04), when we 
adopt their model (see Table 1 of P04). The results are plotted in 
Figure[8] We did not find any significant differences between the at- 
tenuation curves in P04 and ours when using the same dust model, 
suggesting that the presence of spiral arms or any structure of sim- 
ilar size in our simulations does not affect the attenuation signifi- 
cantly. 

It is important to notice that none of the above models, for 
values of dust-to-star scale length/height ratios and central optical 
depths normally observed in spiral galaxies, show attenuations that 
vary linearly as a function of log(a/6). Ferrara et al. ( 1999) found 
that only values of dust-to-star scaleheight ratios of about 2.5 re- 
sulted in an att enuation dependence on incl ination comparable to 
that observed bv lde Vaucouleurs et alj il99lh . but scaleheight ratios 
of this magnitude are far from what observations suggest. Thus, the 



disagreement regarding the behaviour of the attenuation as a func- 
tion of inclination is not limited to our models and it is clear that 
further work is needed to understand this discrepancy. 

In the near future, the images that have been generated with 
StJNRISE for this work will be analysed in order to quantify the 
changes in the apparent scalelengths, central surface brightness, 
and axis ratios due to the presence of dust , in a similar manner 
to what was done bv lMoUenhoff et alj ^20061) . 



6 SUMMARY AND CONCLUSIONS 
6.1 Summary 

Simulations of dust attenuation were performed in eight differ- 
ent SPH simulations of isolated spiral galaxies, assuming different 
metallicity gradients. The three metallicity gradients explored were 
—0.05 dex/kpc, as suggested by Z94, the metallicity gradient re- 
quired to obtain a dust scalelength that is 1 .4 times greater than the 
stellar scalelength, as was observed by X99, and metallicities ob- 
tained from the dependence of the metallicity gradient on luminos- 
ity observed by MA04. We compared the dependence of attenua- 
tion on inclination resulting from these gradients with observations 
by G95, T98, and M03, finding that our simulations do not follow 
a linear dependence on log(a/6). The non-linear fits adopted by 
M03 from their observations of colour dependence on inclination 
were also compared to the results of the simulations, giving a more 
satisfactory match to the shape. The preferred gradients adopted 
from the above comparisons were used to obtain a relation between 
the relative face-on to edge-on attenuation and luminosity, finding 
that it followed a rather linear behaviour that closely matches the 
observations of G95, and T98. Finally, we compare our results to 
previous models of galactic disks where structure such as spiral 
arms are disregarded, showing that the curves obtained from our 
simulations can be explained in terms of the same parameters (no 
parameters to account for large structure resolved by our simula- 
tions are necessary), and that the shape at intermediate inclinations 
is very similar for both the models and our simulations. 



6.2 Conclusions 

1. When the radial metallicity gradient implied by the results of 
X99 is assumed, the observed amplitudes of relative attenuation 
from face-on to edge-on in our most luminous galaxies agree fairly 
well with the results of the simulations. The fainter galaxies of the 
set required less steep gradients, and we adopted the gradients from 
MA04 for the GO, Gl, and G2. This is not surprising considering 
that the X99 results were obtained only for big spiral galaxies, and 
may not be valid for smaller spiral galaxies. 

2. Results for the ratios of attenuations between different 
bands in Table [3] show that there is always some attenuation in the 
near infrared (K-band), with an inclination dependence such that 
the attenuation in the K-band certainly is significant at high incli- 
nations. Our results show that the relative to face-on attenuation in 
the K-band can be as high as 30% of the attenuation in the B-band 
when the Sbc model is seen edge-on, enough to introduce bias of 
about 15% if this is neglected. 

3. Reproducing observed edge-on attenuations and their de- 
pendence on luminosity does not imply that the simulations cor- 
rectly describe the extinction effects across all inclinations. The dis- 
agreement in the shape of the attenuation vs. inclination relations 
between observations and simulations suggests that either a simple 
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Figure 8. Normalised attenuation curves for our Sbc model when the dust model adopted by'Pierini et al. is used in our simulation (left), and for the 

P04 homogeneous face-on model with r = 2 (right). We can see that there are no significant differences that can be related to the structure resolved in our 
hydrodynamic simulations. Figure 18 of P04 shows the dependence of their model with inclination, but only for r = 1, thus we have not plotted these results 
here. 



linear extinction law of the form A\ — -yx log(a/&) is only a rough 
approximation to real galaxies, or that our simulations are not good 
models of observed galaxies. It is difficult to know whether the lin- 
ear relations given by G95 and T98 are the result of uncertainty 
in the observations, or if the relations obtained in the simulations 
really are excluded by the data. The most recent study by M03, 
which also has the largest sample, suggests that a linear extinction 
law is not appropriate, and they propose both a bilinear law and a 
quadratic law. The quadratic fits obtained by M03 have the same 
shape at intermediate inclinations as our simulations. We have too 
few simulated galaxies to perform an adequate statistical fit as has 
been done by M03, hence we cannot adequately compare with the 
amplitude given by these quadratic fits. However, we note that a 
good way to describe the attenuation dependence on inclination 
found in the set of simulations performed here is by a quadratic 
law of the form 

AM = B + C\og{a/b) + D [\og{a/h)f , (14) 

as proposed by [Masters et alj ( |2003|) . Other models obtain a shape 
similar to ours, so this disagreement is not unique to our simula- 
tions 

4. The bolometric attenuations in our simulations as measured 
by the infrared-to-ultraviolet flux ratio are consistent with, but pos- 
sibly on the low side of, the range observed for galaxies of corre- 
sponding morphological type in the SINGS sample by D07. How- 
ever, the bolometric attenuations in the simulation s are significantly 
larger than those obtained in an older study by IPopescu & Tuffs! 

This indicates that the two observational studies are in con- 
flict, regardless of the uncertainty of the simulated results, as it will 
be difficult to match both results simultaneously. 

5. The attenuation vs. inclination relation obtained from mod- 
els that do not take account of spiral structure is very similar in 
shape (a quadratic relation) to our calculations using hydrodynamic 
galaxy models, and there is no evidence to conclude that the in- 
clusion of complex geometries like spiral arms makes a difference 



in the attenuation as a function of inclination nor wavelength. Of 
course, this is not true for highly disturbed systems like merging 
galaxies. 
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